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9^' ABSTRACT 
< 

We obtain self-similar solutions that describe the dynamics of a self- 

^ ■ gravitating, rotating, viscous system. We use simplifying assumptions but ex- 

^ . phcitly include viscosity and the cooling due to the dissipation of energy. By 

O ' assuming that the turbulent dissipation of energy is as power law of the density 



OO 
O 



and the speed Vrms and for a power-law dependence of viscosity on the den- 
sity, pressure, and rotational velocity, we investigate turbulent cooling flows. It 
has been shown that for the cylindrically and the spherically cooling flows the 



O 

P_|. similarity indices are the same, and they depend only on the exponents of the 

Q . dissipation rate and the viscosity model. Depending on the values of the expo- 

nents, which the mechanisms of the dissipation and viscosity determine them, we 
. may have solutions with different general physical properties. The conservation 

of the total mass and the angular momentum of the system strongly depends on 
^ , the mechanism of energy dissipation and the viscosity model. 

Subject headings: hydrodynamics - instabilities - stars: formation - turbulence. 



1. Introduction 

Detailed observations of the interstellar medium (ISM) or even the intergalactic medium 
have highlighted the need to provide a description that accounts for the turbulent pressure, 
thermal pressure, magnetic fields, rotation, and even stars themselves. A universal theory, 
describing the complex structures of the ISM, is far from complete and remains a challenge 
for the future. However, by dividing the ISM according to its properties, it is possible to 
present satisfactory theories for particular types of ISM. Gravitation, cooling, turbulence, 
and magnetic fields produce variations in the properties of the ISM such as the density and 
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the temperature which in turn define the structure of the ISM. Such structures include dense 
cores (e.g., Lee, Myers, & Tafalla 1999), filamentary clouds (e.g., Schneider & Elmegreen 
1979; Harjunpaa et al. 1999), and even disks (e.g., Padgett et al. 1999). Impressive theoret- 
ical progress has been made on the properties and evolution of these structures during last 
years under simplifying assumptions (e.g., Li 1998; Fiege & Pudritz 2000; Tsuribe 1999). 
Most of the theoretical models of interstellar clouds and clumps assume static or stationary 
configurations, in which there is an equilibrium between the self- gravity, centrifugal force 
and some forms of internal energy in the cloud (e.g., Bertoldi & McKee 1992; Chieze 1987; 
Vazquez-Semadeni & Gazol 1995; GaUi et al. 2001; Tomisaka, Ikeuchi, & Nakamura 1988; 
Shadmehri & Ghanbari 2001a). 

Shu, Adams, & Lizano (1987) proposed a four stage scenario for the formation of an 
isolated low-mass star: (a) quasi-static formation and evolution of a molecular cloud core 
by ambipolar diffusion; (b) dynamical collapse of the core to a protostar and circumstellar 
disk; (c) breakout of a powerful bipolar outflow; and (d) clearing of the circumstellar en- 
velope to reveal a pre-main-sequence star. Such theories of isolated star formation assert 
that gravitational collapse occurs onto a thermally supported core and motions are quasi- 
hydrostatic until very late times (e.g., Shu, Adams, & Lizano 1987; Li 1998). Nevertheless, 
recent work has lead to doubts as to whether their initial condition is a reasonable starting 
point. Rather than being quasi-static, the hierarchical and clumpy structures of ISM have 
been proposed to result from turbulence (e.g., Falgarone, PhiUips, & Walker 1991; Elmegreen 
1999). Some authors examine the idea that ISM may not be best visuahzed as a system of 
discrete clouds and conclude that the structures in the ISM may form as density fluctuations 
induced by large-scale interstellar turbulence (e.g., Ballesteros-Paredes, Vazquez-Semadeni, 
& Scalo 1999). However, observations of infall at small scales in isolated cores show that 
quasi-static course do in fact exist (Williams et al. 1999), although it seems that for ex- 
plaining large-scale motions we need an alternative explanation (e.g., Myers & Lazarian 
1998). 

Recently, we presented a non- Jeans scenario for star formation in self-gravitating fila- 
mentary clouds (Shadmehri & Ghanbari 2001b; hereafter SG), extending the work of Meer- 
son, Megged, & Tajima (1996) on spherical clouds. We studied quasi- hydrostatic cooling 
fiows in these filamentary clouds. We parameterized the cooling function as a power-law in 
temperature (A oc p^T^), and showed that the filament experiences radiative condensation. 
Furthermore, since the exponent in temperature is a free parameter, this cooling function 
can also qualitatively represents turbulent energy dissipation. That is we can use this cool- 
ing function to represent the dissipation rate of turbulent energy by identifying the velocity 
dispersion with the temperature (see, e.g., Mac Low 1999). Thus, our results also apply to 
a turbulent filament, although we did not investigate this explicitly. 
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Myers & Lazarian (1998) have suggested that inward, subsonic flows arise from the local 
dissipation of turbulent motions in molecular clouds. These turbulent cooling flows result 
from localized dissipation of turbulence in molecular clouds, as radiative cooling flows result 
from radiative cooling in clusters of galaxies (e.g., Cowie & Binney 1977). Myers & Lazarian 
(1998) considered only one type of turbulence (a linear superposition of MHD waves) and 
only one dissipation mechanism (ion- neutral friction). They also ignored the role of self- 
gravity. On the other hand, there are competing explanations for the origin of the clumps 
in the molecular clouds. A more accepted explanation is that the clumps continuously form 
and disperse simply because the gas is supersonically turbulent (e.g., Scalo 1990; Falgarone 
& Phillips 1990). Recent numerical studies of turbulence have brought a new understanding 
of the physics of this complex phenomenon. For example, recently Padoan et al. (2001) 
suggested that because of the turbulent nature of supersonic motions in molecular clouds, 
the dense structures such as filaments and clumps are formed by shocks in a turbulent flow. 

In this paper, we extend the work of SG to study turbulent cooling flows in self- 
gravitating, rotating, viscous clouds. We adopt an analytic approach and do not perform 
extensive numerical computations. Viscosity has an important role in the behavior of the 
cloud. Of course, molecular viscosity is totally negligible in the ISM compared to turbulent 
viscosity. This presents a problem as no comprehensive treatment of turbulent viscosity ex- 
ists. For simphcity, we shall use the well-known a— prescription (Shakura & Sunyaev 1973) 
which has been used for modelling accretion disks and also a more general simplified model 
for viscosity. This study is a first step (at least, qualitatively) toward understanding the 
role of turbulent coohng flows in self-gravitating clouds. Clearly, turbulence is a complex 
phenomenon in the ISM and the results of this investigation arc just approximations to the 
behavior of real turbulent cooling flows. In general, our treatment will be incomplete because 
we omit turbulent eddies with scales larger than the size of the cloud. Thus, we consider 
microturbulcncc, and this can only provide a flrst approximation to the real dynamics of 
ISM clouds. Nevertheless, our simple model provides useful information on the importance 
of turbulent cooling flows in star forming regions. 



2. General Formulation 

We shall seek similarity solutions of the set of equations describing turbulent cooling 
flows in sclf-gravitating clouds. The simplest method of dealing with turbulence involves 
the assumption that the Navier-Stokes equation holds on each physical scale, with turbu- 
lence manifesting itself via renormalized viscosity and heat conduction coefficients that are 
generally not constant. We assume spherical (or cylindrical) symmetry and neglect heat con- 
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duction. For constructing the model, the continuity equation, the momentum equation in 
the radial direction, the equation of angular momentum conservation, the Poisson equation, 
and the energy equation can be re-written explicitly in the following way: 

dp , 1 d 



dt r^~^ dr 



(r^-'pv,) = 0, 



dvr dvr I dp 9* _ vl 
dt ^ dr pdr dr r ' 

dt ^""^ dr ~ r^-^pdr^"^^' dr>' 

where p, w,., v^, p, and ^ denote the gas density, radial velocity, rotational velocity, pres- 
sure, and gravitational potential, respectively. Also, Q = ^ is the angular velocity and A 
determines the dimensionality of the model, with A being 2 in cylindrical geometry and 3 in 
spherical geometry. 

To close this system of equations, the form of and viscosity (u) cooling function (A) must 
be known. As was mentioned above, our understanding of turbulent viscosity is incomplete, 
and for this reason we adopt an empirical prescription. For example, the a— prescription 
assumes that the molecular viscosity can be replaced by an isotropic turbulent viscosity, and 
so 

where a is a constant (Shakura & Sunyaev 1973). It has been shown that a is not in general 
constant, but rather depends on a number of factors (see, e.g., Brandenburg 1998). However, 
for simplicity, we will assume that the parameter a is constant. Also, we consider a more 
general viscosity model than a— model, viz., 

u = ap^pf^n", (7) 

where r, /?, and rj are three arbitrary exponents. A viscosity model similar to this form 
has been used by Begelman & Meier (1982). By imposing additional requirements like 
conservation of the total mass or angular momentum of the system, we shall find relations 
between these three exponents. It is clear that the a— model corresponds to 



T = l,/3 = -l,?7 = -l, 



(8) 
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and is one member of this family. Note that the constant a in equation (7) is in general 
different from the a in equation (6); they are the same only when the exponents have values 
given by equation (8). We use a as a free parameter to study the effect of viscosity 

For dissipation, we could use the Mac Low (1999) result locally, which would give a 
dissipation time that would be always comparable to the local crossing time. In reality, 
the dissipation rate is not just a function of temperature. The T varies all over because of 
shocks, and most of the dissipation is in shocks , not in a uniform thermal medium. So, if 
we have some characteristic scale, like the radius around a filament or the distance from the 
center of a sheet, we would take pv'^.^^^jL as the dissipation rate for speed Vrms- Then L = r 
of where we consider in the filament, for example. Also, some authors have found power-law 
relations between Vrms ^-nd L by doing numerical simulations (e.g., Vazquez-Semadeni & 
Avila- Reese 2001). We note that it is possible by knowing the Vrms, define the temperature 
T. Neglecting any viscous and external heating, however, it is possible to fit the cooling 
function (dissipation rate) by a power law, 

A(p, T) = /l,/T^ (9) 

where and e are constants. This general form enable us to study the behavior of the system 
by changing the values of e and also A^. In fact, the mechanism of turbulent dissipation 
determines the values of these parameters. 

In certain optically thin systems, we can also use equation (9) to approximate the 
radiative cooling function. Radiative cooling almost always involves a two body process, 
and as such also depends of the square of the density. The exponent of temperature e 
depends on the regime under consideration. In molecular clouds, cooling is dominated by 
dust or CO line emission, and the range of e is from 1.5 to 3 (Goldsmith & Langer 1978). 
Also, for hot plasmas of cosmic composition, we have e — | for T > lO^i^ when free-free 
emission dominates, and e ii: — | for lO^i^ < T < 10^ K when fine cooling dominates (Gaets $ 
Salpeter 1983). This approach is useful for quahtative investigation of the effects of radiative 
cooling, but we must be aware of its limitations. 

We are interested in quasi-hydrostatic flows in the clouds. For quasi-hydrostatic flows, 
equation (2) becomes 

if! + §^^4 (10) 

par or r 

To simplify the equations, we make the following substitutions: 

p ^ pp,p ^ pp,^ T TT, Vr,if, vvr,y,, r ^fr,t ^ it, (11) 



where 
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Under these transformations, equations (1) and (10) do not change, but the other equations 
become 

+ = -^^(pV*'n"r'»f). (13) 

at or mr^'^por or 

P^a^^' ^^^^ 

Notice the equations are not invariant under the transformation t ^ —t, Vr ^ —Vr, and 
— > —Vy;, so the solutions can not be time-reversible. Thus, these equations describe a 
nonlinear evolving system. 



3. Similcirity Solutions 
3.1. Analysis 

In following the nonlinear evolution of dynamically evolving systems, the technique of 
self-similar analysis is useful, as it allows a set of partial differential equations, such as those 
above, to be transformed into a set of ordinary differential equations. A similarity solution, 
although constituting only a limited part of the problem, is often useful in understanding the 
basic behaviour of the system. In order to seek similarity solutions to the above equations, 
we introduce a similarity variable ^ as 

(16) 



and assume that each physical quantity is given by the following form: 

(17) 

(18) 
(19) 
(20) 
(21) 



p(r, t) : 


= {to- 


-trR{0, 


p(r, t) ■■ 


^{to- 


-trp{0, 


Vr{r,t) 


= {to 


-trv{o, 


v^{r,t) 


= {to 


-trm, 


*(r, t) 


= {to 


-trs{o, 
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where the exponents n, ui, U2, 1^3, 1^4, and are constants which must be determined and 
t < to- By substituting the above equations into equations (1), (10), (13), (14), and (15), we 
obtain the general results: 

2n(l - e) - 1 2n(2 - e) - 2 ^ 2n-l 2n-l 

i/i = ,1/2 = ,i/3 = n - 1,1/4 = — ,i^5 = , (22) 

e e 2e e 

and, 

+ + (23) 
(?7 + 2/3 + 2T + 2)e- (7^ + 2/? + 4t)' ^ ^ 

In the specific case of the ct— model, we have t — 1, (3 — rj — —1, which gives: 

2e-l . 3-2e 1 1 1 

^=2(731)'-! = = ^,-3 = ^(731)'-^ = ^(731)'-^ = 731- (24) 

It is interesting that these similarity indices are independent of the dimensionality of the 
system, i.e. A, but depend strongly on the exponents of the viscosity and the turbulent 
cooling function. Thus, in both cylindrical and spherical just the mechanisms of viscosity 
and dissipation of energy are important. 

We can write the equations for the dependence of the physical quantities on the similarity 
variable as: 

-''.« + '^f + |^|«^-'«V) = 0. (25) 

1 dP dS $2 

- •'2P + + '''^ + + ri'"P' = 0. (28) 

This is a system of non-linear ordinary differential equations which has some critical 
points. Since the similarity indices arc known, it is possible to study the general physical 
properties of the solutions. Wc sec that the solution for each physical quantity retains a 
similar form as the flow evolves, but the characteristic length-scale of the flow increases 
proportionally to (to — t)"', where equation (23) gives the value of n. As we require that n 
be finite and nonzero, this relation shows that e may have any value except e = ^^2^^tr+2 
or € — P + 2t + The sign of n determines whether the cloud expands or collapses; if n is 
positive the radius of the cloud Tc decreases with time and the central density Pc increases 
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with time. In the a— model, n depends only on e, and is negative for e between | and 1, 
and positive otherwise. Also, in this case, the density at the center of the cloud increases in 
proportion to {to — irrespective of the value of e. 

We now consider boundary conditions, such as conservation of the total mass, angular 
momentum and constant external pressure. Conservation of mass is an important physical 
constraint, and, as we shall show, fortuitously simplifies the equations and allows us to 
obtain some analytical results. The integrals representing the total mass (M) and angular 
momentum (J) are the following: 

M= fyr^-^dr, (30) 



J 



= J pUr^+^dr. (31) 



Thus, the total mass M (nondimensional) is proportional to {to — t)"'^'^^'", and the total 
angular momentum J (nondimensional) proportion to {to — tY'^^"^'^^^'^^^'^ . As these relations 
show the conservation of M or J depends on the dimensionality of the system and the 
exponents of viscosity and cooling rate. First, we study the properties of the solutions for 
the a— model. There is no value of e for which the mass of a cylindrical cloud to be conserved 
and for spherical clouds the mass of the cloud is conserved only for e = — |- The total angular 
momentum is conserved with e = 2{x'-i) ■ However, when the total mass of a spherical cloud 
is conserved (e = — |), the total angular momentum decreases: J cc {to — t)^. Another 
interesting case is a cylindrical cloud with general form of viscosity, in which the constraint 
of constant total mass gives 

77 + 2/3 + 2t = 0. (32) 

Note that this constraint involves only the indices of the viscosity model and not the in- 
dex of the dissipation model, e. Thus, this constraint is independent of the mechanism of 
dissipation. Under this constraint, equations (22) give the similarity indices: 

1 1 

n = -, i/i = -1, i/2 = -1, = = 0, i/5 = 0. (33) 

Also, in this case, the total angular momentum decreases: J a {to — t)^. 

Constant external pressure is another interesting case in which the cloud is surrounded 
by an ambient gas: p{r — > oo) = const. Of course, the total mass need not be conserved, as 
inflow or outflow of the gas is allowed. In the a— model this corresponds to e = |, regardless 
of dimensionality, and the total mass and the angular momentum of the system decreases: 
M a (to - t)'(^-'^ and J a (to - t)^^^+^\ 
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Prom the above, it is clear that the solutions show different behaviors depending on 
the values of the exponents of the dissipation of energy and the viscosity (that is, on the 
dominating mechanisms of dissipation and viscosity). Also, these results are independent of 
the magnitude of the turbulent viscosity. 



3.2. Numerical solutions 

Once e, r), (3, and r are selected, we can solve the set of ordinary differential equations 
(25)-(29). The full range of possibilities becomes enormous if we regard all the exponents 
of the viscosity model and the cooling rate as free parameters. For this reason, we shall 
investigate only the case of constant total mass. First, we study cylindrical cooling flows 
with constant total mass, in which the constraint given by equation (32) must be satisfied. 
If we require a finite density and zero velocity at the cloud center, equation (25) gives 

V = -nC, (34) 

where n = |. Substituting this equation into equation (28), gives an algebraic relation 
between the similarity functions of the pressure and the density: 

P = PoR", (35) 

where Pq = (7 — 1)~ and q = It is interesting that this relation for turbulent and 

rotating cooling flows is the same as the one derived by SG for radiative cooling flows. Note 
that there is no critical value for 7, because for all values of 7 > 1, the relation between the 
density and the pressure is well-defined. Also, e may have any values except e — 1 and 2. 
Going back to equations (27) and (29) and using this result, we have: 

|«3-,«...«*,|,£))^_^5.«.0, (36) 

|(4^-5P„eif-^^)-eii=0. (37) 
We have to solve these equations subject to the boundary conditions: 

R{C = 0) = 1, (^),=o = 0, $(e = 0) = 0, (^),=o = 0. (38) 



Of course, there is another constraint: conservation of mass. Since we are not interested in 
the exact value of the mass of the system, we don't parameterize the solutions by the total 
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mass. In fact, we integrate the equations by the Runge-Kutta method to seek solutions sat- 
isfying the boundary conditions and then we can calculate the mass of the system. However, 
it could be simply checked that ii q = then the above equations are invariant under a 
homology transformation. Indeed, if R{C) ^ solutions of the equations, 

1 2g-3 

then Ai-'^R{A^) and A'-'^i-'-'') ^(A^) arc also a set of solutions, where A is an arbitrary con- 
stant. Using this transformation, it is possible to find A so that the solutions are satisfying 
the boundary conditions and the constraint of constant mass. Thus, if r and q are satisfying 
the relation q = we can parameterize the solutions by the line mass of the cylindrical 
cloud. 

In Figure 1, we represent the distributions of the density -R(C) in the self-similar space 
for 7 = |, m = 4.0, a = 0.007, r = |, 77 = —1, /3 = and different values of e = 3, 4, 
—3, —4. They are obtained by solving the equations with the fourth-order Runge-Kutta 
method. This Figure shows that the general behaviors of the similarity density function 
R{^) for different values of e are almost the same. In fact, all the solutions tend to zero as 
^ — > 00. It means that we can define the radius of the cloud by the condition R{^c) — 0, 
where Tc = ^c{to —t)^ (see equation (33)). 

Figures 2 and 3 show the behavior of the similarity rotational velocity function $(^) for 
the parameters the same as Figure 1. From Figure 1 and 2 we see that whole of the cloud 
is rotating for e = 3 and 4. But Figure 3 shows that only part of the cloud is rotating for 
e = —3 and —4. Also, each profile shows a maximum value for the rotational velocity. 

For the a— model viscosity, we showed that only for e = — | the spherically cooling 
flows describe a system with constant total mass and this value of e corresponds to free-free 
emission for cosmic compositions. Figure 4 shows the density profile of spherically symmetric 
fiows for different values of a. As the Figure shows the profile of R{C,) hardly depends on 
the value of a. It means that the distribution of the mass in spherically symmetric cooling 
flows is independent of the magnitude of viscosity. But as Figure 5 shows, the distribution 
of angular momentum in such systems strongly depends on the magnitude of viscosity. In 
Figure 5, we can see the normalized rotational velocity $({) for different values of a and 
it shows as the magnitude of viscosity increases, the maximum value of the profile of $(^) 
increases. 



4. Discussion and Conclusions 

We have explored in this paper a new set of similarity solutions of the equations relevant 
to cooling flows of self-gravitating, rotating, viscous systems. Analogously to the case with 
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no rotation and viscosity (SG; Meerson, Megged, & Tajima 1996), we were able to obtain 
similarity solutions which describe cooling flows in rotating and viscous systems. However, 
we did not investigate all of the solutions of the system. In fact, wc restrict ourselves to 
the quasi-hydrostatic flows with constant total mass under a wide range of the viscosity and 
the cooling models. When the central density is high and the velocities small, we can apply 
our solutions which are regular at the center. Even if we solve the equations without the 
quasi-steady assumption by including the time dependent terms in the momentum equation, 
we shall find the same similarity indices as equation (22). 

It may seems that physical constraints on our solutions, such as mass or angular mo- 
mentum conservation, limit the acceptable range of similarity indices. For example, the 
a— viscosity law is incompatible with mass conservation in our models. But it does not mean 
that the self-similarity assumption is not a good assumption. As for the viscosity model, 
we mentioned that since within the turbulent viscosity assumption the main difficulty is to 
point out a plausible source of turbulence, empirical viscosity prescriptions (e.g., a— model) 
have been broadly used. Although the a— model enables a reasonable global description of 
stationary, thin and non self-gravitating disks (Frank et al., 1992), we must note that the 
a— model is only an empirical prescription. Few improvements of the a— prescription have 
been made (Narayan 1992; Narayan et al, 1994; Godon 1995) and it has been shown that 
in general a is not constant (Brandenburg 1998). Thus, just according to considerations 
about viscosity models, we can not investigate the validity of our solutions. Due to these 
facts, we introduced a general form for viscosity law, i.e. equation (7). However, one should 
demonstrate that the acceptable viscosity laws have reasonable dependencies on physical 
quantities. But such study is out of the scope of this paper. 

Recently, Elmegreen (1999) investigated the role of wave-driven turbulence in interstellar 
clouds. In his study, a dense region forms because of a convergence and high pressure from 
external magnetic waves and the thermal cooling that follows at the compressed interface. 
However, turbulent dissipation at the interface is not as important as thermal cooling for this 
density enhancement (Elmegreen 1999). Thus, if in our model we consider the dissipation 
rate simply as thermal cooling function (not as turbulent dissipation rate), probably it will 
be possible to apply the results of this study for the compressed interfaces in Elmegreen's 
model. 

In fact, some cosmological structures, diffuse HI mediums, and molecular clouds, all 
of which show some cooling properties. Thus, our solutions might have applications to the 
formation of these structures at early stages. In this regard, this work is an extension of the 
same cooling flow problem considered in SG and Meerson, Megged, & Tajima (1996) in order 
to study the effects of viscosity and rotation on cooling ffows. For filamentary clouds, SG 
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showed that if the mass of the system is conserved, then the density at the center increases 
{pc oc (to — t)~^) and the radius of the filament decreases oc (to ~ ^)^), irrespective 
of the e, that is, the mechanism of cooling. Wc showed that if we consider the rotation 
and viscosity, conservation of the total mass of the system depends on the mechanism of 
cooling and viscosity model. For example for the a— model, all the cooling mechanisms 
lead to a system with nonconstant mass. But for the other models of viscosity which the 
exponents satisfy into equation (32), the total mass of the system is conserved and we have: 
Pc cc (to — t)~^ and Vc oc [Iq — t)^. Also, for spherical coohng flow the value of e = — | 
corresponds to constant mass. In this case, we showed that pc oc (to — 1)~^ and Tc oc (to — t) 3 . 
However, Meerson, Megged, & Tajima (1996) showed that for nonviscous and nonrotating 
spherical cooling flows, the behaviors of and pc depend on e and the mass of system is 

1 3 

conserved: Tc oc (to — t)^ and pc oc (to — t)~2+;. 

A dispersion relation has been derived for gravitational instabilities in a medium with 
cloud collisional cooling by Elmegreen (1989). The cooling function due to the cloud-cloud 
collisions with an isotropic, Maxwellian distribution of cloud velocities is a power-law of the 
density and the temperature (Elmegreen 1987), similar to the form of our cooling rate in 
this paper. Thus, in a cloudy medium without much star formation activity (i.e., without 
heating), the energy dissipation is given by equation (9). Elmegreen showed that the regions 
with such conditions will clump into cloud complexes on a variety of scales. Since on the 
largest scales, the rotation and shear of the Galaxy become important, our results describe 
the general properties of structure formation in such regions. 

This simple model has nevertheless some limitations. We have not included explicit 
heating terms in our energy equation, in order to minimize the number of parameters. How- 
ever, we can consider additional terms in the governing equations. For example, one could 
invoke the flow inertia in the momentum equation, or heat conduction in the energy equa- 
tion. Although the similarity indices would be selected in these cases, we can consider only 
limited values of the exponents of the cooling function and the viscosity model. Also, the 
assumption of similarity solutions, could be dropped, although we should not forget the com- 
plex nature of the governing partial differential equations. Also, we assumed a power law for 
the form of cooling rate and viscosity and the results depend strongly on the form of these 
functions. However, as we mentioned in SG, it seems that we can use piece-wise solutions. It 
means that for each part of the evolution of the solutions, we can flnd the suitable exponents 
and then join one set of solutions to the another one. Thus, it will be possible to consider 
a much wider range of cooling functions and viscosity models with the same analytical-type 
of solution. 

This work was completed while one of the authors (MS) was visiting Instituto de 
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Fig. 1. — Similarity solutions for the cylindrically symmetric quasi- hydrostatic flows for 
7 = |, m = 4.0, a = 0.007, r = |, = —1, and P = 0. Shown are the normalized density 
R{C) for different values of e = —4 (solid curve), —3 (dotted curve), 4 (dashed curve ), and 
3 (long- dashed curve). All the solutions tend to zero as ^ — > oo. 
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Fig. 2. — Similarity solutions for the cylindrically symmetric quasi-hydrostatic flows for 
7 = |, m = 4.0, a — 0.007, t — \^r\ — —1, and /3 = 0. Shown are the normahzed rotational 
velocity $(^) for different values of e = 4 {solid curve), 3 {dotted curve). 
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Fig. 3. — Same as Figure 2, but for different values of e = 



—3 (dashed) and —4 (long- dashed). 
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Fig. 4. — Similarity solutions for the spherically symmetric flows for 7 = |, m = 4.0, 
e = — |, and a— model viscosity. Shown are the normalized density -R(0 different values 
of q; = 0.005 {solid curve), 0.05 {dotted curve), and 0.1 {dashed curve). We can see that R{C) 
hardly depends on a. 
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Fig. 5. — Similarity solutions for the spherically symmetric flows for 7 = |, m — 4.0, 
e — —J, and a— model viscosity. Shown are the normalized rotational velocity $(^) for 
different values of a = 0.005 {solid curve), 0.003 {dotted curve), and 0.001 {dashed curve). 



